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Abstract: The overshooting convective motions in the solar photosphere are fre- 
quently proposed as the source for the excitation of Alfven waves. However, the pho- 
tosphere is a) very weakly ionized, and, b) the dynamics of the plasma particles in this 
region is heavily influenced by the plasma-neutral collisions. The purpose of this work is 
to check the consequences of these two facts on the above scenario and their effects on 
the electromagnetic waves. It is shown that the ions and electrons in the photosphere are 
both un-magnetized; their collision frequency with neutrals is much larger than the gyro- 
frequency. This implies that eventual Alfven-type electromagnetic perturbations must 
involve the neutrals as well. This has the following serious consequences: i) in the pres- 
ence of perturbations, the whole fluid (plasma -|- neutrals) moves; ii) the Alfven velocity 
includes the total (plasma -|- neutrals) density and is thus considerably smaller compared 
to the collision-less case; iii) the perturbed velocity of a unit volume, which now includes 
both plasma and neutrals, becomes much smaller compared to the ideal (collision-less) 
case; and iv) the corresponding wave energy flux for the given parameters becomes much 
smaller compared to the ideal case. 
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1 Introduction 



In a weakly ionized but highly collisional medium, a propagating Alfven wave also involves 
the motion of the neutrals that are present in the medium. This is due to the friction 
between charged particles and neutrals. The effect has been described in the hterature, 
first by Tanenbaum and Mintzer (1962), Woods (1962), Jephcott and Stocker (1962), and 
in many subsequent works, e.g., Kulsrud and Pierce (1969), Pudritz (1990), Haerendel 
(1992), De Pontieu and Haerendel (1998), Watts and Hanna (2004). This fundamental 
result is valid for any weakly ionized plasma, including the plasma in the lower solar 
atmosphere. 

The Alfven wave has been a very popular tool in the scenarios and models dealing 
with the heating of upper solar atmosphere. A necessary ingredient in such models is an 
efficient and abundant source for the excitation of these waves, which acts permanently 
and generates waves throughout the solar atmosphere. Very frequently it is assumed that 
the omnipresent overshooting convective motions in the photosphere could serve for this 
purpose. The amount of thermal energy per unit volume in the solar corona is in fact 
extraordinarily small in comparison with the lower (and much colder) layers of the solar 
atmosphere. This is due to the rapidly decreasing density with altitude. On the other 
hand, the complete photosphere is covered by overshooting convective gas motions with 
typical velocities of about 0.5 km/s, that may go up to 2 km/s. The kinetic energy 
per cubic meter stored in this macroscopic motion of a mainly neutral gas exceeds for 
several orders of magnitude the internal energy density in the corona. Clearly, only a tiny 
fraction of the convective kinetic energy of the neutral gas would be sufficient to heat 
the higher layers to the given temperatures. Such a scenario is attractive in view of the 
fact that this macroscopic motion in the lower atmosphere is permanent and widespread 
throughout the solar surface. However, the photosphere is very weakly ionized and it is 
also a strongly collisional mixture of the tiny plasma component and the predominantly 
neutral (uncharged) gas. 

The energy flux of the Alfven waves is given by minivfca/2^ where is the Alfven 
velocity and Vi is the perturbed velocity of ions involved in the oscillations. Typically, in 
the estimate of the flux in the photosphere, this perturbed velocity taken is of the same 
order as the macroscopic convective motion mentioned above (HoUweg 1981). 

In the present work, we focus on the physics involved in the propagation of the Alfven 
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Table 1: Collision frequencies (in Hz) and magnetization ratio of electrons and protons 
in the photosphere for two altitudes h (in km) and for the magnetic field Bq — 10"^ T. 



h 

















1.6 ■ 10^ 


5-10^ 


1.3 ■ 10^° 


1.5 ■ 10^ 


6 ■ 10-4 


1.1 • 10-1 


250 


2.6 • 10^ 


3.8 • 10^ 


2.2 • 10^ 


1.2 • 10^ 


3.6 • 10-3 


7.3-10-1 



wave in a weakly ionized plasma like the solar photosphere. Using simple and reliable 
physical arguments and widely accepted plasma theory, we discuss the flux of the Alfven 
waves under these circumstances. It will be shown that, if we assume the existence of the 
necessary electromagnetic perturbations in such a weakly ionized medium, the energy flux 
of the waves is in fact much lower compared to what is usually expected from estimates 
based on ideal magnetohydrodynamics. This is due to the fact that the photospheric 
gas dynamics is heavily influenced by collisions. More precisely, in the presence of some 
accidental electromagnetic perturbations, which in the first step involve plasma species 
(electrons and ions) only, the neutral atoms respond to these electromagnetic perturba- 
tions due to the strong friction. This, and the fact that the ionization ratio is rather small 
(viz. of the order of 10-^), results in very small amplitudes of the perturbed velocity of 
the total plasma-gas fluid. 

2 Physics of weakly ionized plasmas 

We introduce here the collision frequencies between charged and uncharged particles Uj^ — 
nnoCjnVrj for j — 6, i, and the formulas (Mitchner and Kruger 1972; Spitzer 1963) for the 
Coulomb coUisions between charged plasma particles: 

+ Vei^2vei= [Aueoi^i: / mef'ieei/ {AT:So)f Lei/ [2,{KTef'^]] , 

uu = [4n,o{7r/m^y/^[e^J{Aneo)]^Lu/[3{KT,f^] . 

All frequencies will be given in Hz and Lei — log[127r£o(£o/?T'io)^^^(/«T'e)^'^^/(eef)] describes 
the Coulomb logarithm. 

Several comments are noteworthy before continuing the derivation. Using the full 
quantum theory as well as the semi-classical approach, the elastic proton-hydrogen [H^ -\- 
H) coUision cross section Uin is calculated by Krstic and Schultz (1999), and its integral 
value at 0.5 eV is about 1.8- 10~i^ m^ for the elastic scattering, and about 10~i^ m^ for the 
momentum transfer. As for the electron-hydrogen {e~ + H) collisions, the collision cross 



section a^n is also temperature dependent and the corresponding values can be found in 
the works of Bedersen and Kieffer (1971), and Zecca et al. (1996). At energies of 0.5 eV 
it is about 3.5 • 10~^^ m^, so that for the elastic scattering we have Oin/cFen — 6. 

On the other hand, here we do not include the inelastic collisions that take place in a 
partially ionized plasma, like in the photosphere. It can be shown (Vranjes and Poedts 

2006) that, in the photosphere, all ions in a unit volume are recombined many times per 
second. The thrcc-body recombination (the process of the type + e~ + e~ ^ H + e~) 
is dominant in this region. At the altitude of h = 500 km, the radiative recombination 
(the process described by + e~ —>■ H + hu) and the three-body recombination are of 
the same order. At higher altitudes, the radiative recombination becomes the leading loss 
effect. At h — 1000 km, it is by a factor 100 larger than the three-body recombination. 

In addition, the charge exchange between the ionized and neutral hydrogen is frequent. 
The cross section (Krstic and Schultz 1999) for the proton-hydrogen charge exchange agx 
at the above given temperatures is about 5.6 • 10^^^ m^, i.e., for hydrogen it is a large 
fraction (~ 0.3) of the realistic elastic scattering cross section cTj^ given above. Note, 
however, that for some other gases, like He, Ne, and Ar, we have a^.^ > (Raizer 
1991), i.e., the charge exchange cross section exceeds the one for the elastic scattering. 
Consequently, due to the inelastic collisions and the charge exchange, neutrals/ions in the 
plasma spend a part of their time in the ionized/neutral state, respectively. As a result, 
the effective collision frequencies are expected to be even higher than the values that we 
shall use. 

Using the data for a quiet Sun model (Vernazza et al. 1981), in Table 1 we summarize 
the values for the electron and proton elastic scattering collision frequencies at two alti- 
tudes (viz. h = km, and h = 250 km) in the solar photosphere (see also Vranjes et al. 

2007) . Here, we have taken Bq = 10"^ T, the corresponding temperatures arc respectively 
T = 6420 K and T = 4780 K, the electron number densities are rio = 6.4 • 10^*^ m-3 and 
no — 2.7 ■ 10^^ m^'^, and the atomic hydrogen number densities are n„o = 1-17 ■ 10^^ m~^ 
and Uno — 2.3 • 10^^ m~^. We assume that the proton and electron number densities are 
equal. It is seen that both protons and electrons are un-magnetized. Note that in Table 1 
the collision frequencies between the plasma species and neutrals are dominant for both 
electrons and ions, compared to the frequencies for Coulomb collisions between charged 
particles. 

It is believed (Priest 1987; Sen and White 1972) that, due to the low temperature. 
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the ions in the lower photosphere are in fact mainly metal ions. Sen & White (1972) have 
assumed that the mean mass of these metal ions is 35 a.u. In that case, due to the rather 
different masses of (metal) ions and neutral (hydrogen) atoms, in calculating the collision 
frequency it is appropriate to use a more accurate formula — nnQ(7mn[^n/ {i^m + 
mn)][9>KT„J {jiji)]^/'^ , where the index m denotes the metal ion, n denotes the neutrals 
(hydrogen), and ^ = 171^171^/ {mm + m„) is the reduced mass. The calculations may be 
inaccurate because the collision cross section (Tmn is not known. As a guess, we take it 
as the value for protons multiplied by mm/nip. Taking the layer h = 250 km, we find 
^mm — 6.4 • 10^ Hz, Vmn — 4 • 10^ Hz, and Vtm — 2.7 • 10^ Hz. Comparing to protons 
from Table 1, the metal ions appear to be even less magnetized, i.e., fl^/u^ — 6.6 • 10~^, 
where i/^ — v^^n + ^mn- At /i = km, we have v^^n — 1-2 • 10^ Hz, i/^„ = 2 • 10^ Hz, 
and D.yn/i'm — 1-3 • 10"^. The mentioned uncertainty in determining (7^„ will clearly not 
substantially change the fact that the ions are un-magnetized. 

3 Physical picture of Alfven waves in a weakly ion- 
ized plasma 

Following standard textbooks (e.g. Chen 1988), in the case of the shear Alfven wave with 
Bq = BqCz, both ion and electron fluids oscillate in the direction of the perturbed magnetic 
field vector Bi = BiCy. This is due to the Ei x Bq drift, which separates neither charges 
nor masses, and the direction of the electric field is determined by the Faraday law. The 
wave is in fact sustained by the additional polarization drift Vpj — {ruj / qjBl)dEi/ dt 

— * 

and the consequent Lorentz force j^ex x Bq, which is again in the y-direction and has a 
proper phase shift. Note that the polarization drift appears as a higher order term due 
to \d/dt\ <^ Jlj. It introduces the ion inertia effects and if it is neglected, then the Alfven 

— * — * 

wave vanishes. The E x B term essentially describes the magnetic field frozen-in property 
of the plasma. The mode is fully described by the wave equation 



VxVxEi^ —E, + ^ji. 
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the momentum equations for ions and electrons 
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-mnnnVni{Vn - Vi) - m„nniyne{Vn - Ve)- (4) 



- UleneiyeniVe " Vn) - meUeiyeiiVe " Vi), (3) 

and the corresponding equation for neutrals 
dt 

Usually, the viscosity may be omitted, and this is valid even for the neutrals. We note 
that the dynamic viscosity coefficient for the atomic hydrogen fi for the layers h = and 
250 km (with the temperatures Tq = 6420 K and To = 4780 K), can be taken (Vargaftik 
et al. 1996) as 6.5 ■ 10~^ Ns/m^ and 5.2 ■ 10~^ Ns/m^, respectively. The corresponding 
kinematic viscosity coefficient rj — ^/{mnnn) is 1.34 m^/s. The ratio rjk^/uni for the 
wave-lengths of interest is very small and the viscosity effects appear negligible in spite 
of such a low ionization. Further, using the momentum conservation in the friction force 
terms yields i/ni — rriinQVin/ {mnnno), and &X, h — 250 km and for = m„ we obtain 
Vni — ?>- 10^ Hz. The thermal effects may also be omitted, as will be explained below. 

In the absence of collisions, the response of a plasma to the magnetic and electric 
field perturbations is instantaneous, and a volume element of the plasma moves in the 
previously described manner. In such an ideal case, the energy flux of the Alfven wave is 
given by 

Fid = minovfca/2. (5) 

— * — * 

Here, Vi is the leading order E x B perturbed ion velocity. Its amplitude is given by 
Vi — Ei/Bq. Using the Faraday law we have Ei — uBi/k, hence Vi — CaBi/Bo. For 
the estimate only, we assume small perturbations of the magnetic field, viz. around 1% 
(a comment on larger perturbations will be given later on). For the parameters ai h = 
250 km, this yields Ca = Bq/ {fj,oniomi) = 1.3 ■ 10^ m/s. Consequently, the perturbed 
plasma (ion) velocity is Vi — 10~^Ca = 1.3 ■ 10^ m/s. The wave energy flux in the ideal 
case, and for rrii — nip, becomes Fid — 5.3 • 10^ J/(m^s). Setting rrii — 35mp yields 
Fid 90 J/(s m2). 

CoUisions may heavily alter the motion of the perturbed electron and ion fluids. The 
plasma response to the electromagnetic (Alfven-type) perturbations in fully and weakly 
ionized plasmas is essentially different from various points of view. Here, we present some 
facts that should help in understanding the physics involved in the description of the 
Alfven waves in the partially ionized plasmas, and in particular in the photosphere. 
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Figure 1: Schematic presentation of the motion of a charged particle in non- magnetized 
plasma. 

1. From Eq. ([2]) we see that the ratio of the Lorentz and the friction forces (for pre- 
dominant ion- neutral collisions and in the case of initially unperturbed neutrals), is 
Qi/uin. For the given photospheric plasma this is ~ 1/10^. 

2. The motion of an un- magnetized charged particle is depicted in Fig. 1. Arrows 
denote the tangential direction at the moment of collision when the particle switches 
to another gyro-orbit with a possibly different velocity (indicated by different gyro- 
radius). A collision occurs after a very tiny fraction (largely exaggerated here) of 
the gyro orbit has been traveled. According to the numbers from the Table 1, the 
particle trajectory along a gyro-orbit around one specific magnetic line is only about 
1/10^ part of the full circle. Hence, in the given case the path of the particle between 
two collisions is nearly a straight line, like in the case when the magnetic field is 
absent. In fact, the particle never makes a full rotation. The motion is similar in the 
fully ionized plasmas, however, there it is related to ion-ion collisions (i.e., viscosity, 
not to friction). 

3. Contrary to the case of fully ionized plasma where all particles in a volume element 
move together due to the given electric field while still colliding with each other, in 
the present case the collisions introduce a limitation. In the beginning only charged 
particles are supposed to move due to the applied electric field, while neutrals have 
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a tendency to stay behind. If ^ ^i,Oi), each plasma particle collides many times 
within the theoretical gyro-rotation, or within the assumed wave oscillation. 

4. Contrary to the viscosity, which is of primary importance for short scale processes, 
the friction is more effective in the opposite limit, i.e., at smaller wave-numbers 
(and also at larger wave-periods) an ion is subject to larger number of collisions 
with neutrals within one oscillation period. 

Below, we discuss the effects of collisions in detail. According to the presently widely 
accepted physical description of Alfven waves in partially ionized plasmas (Kulsrud and 
Pierce 1969; Haerendel 1992; De Pontieu and Haerendel 1998; Pecseli and Engvold 2000), 
the effects of neutrals are the following: a) for a relatively small amount of neutrals (or 
for high frequency short wavelengths), the damping of the mode is proportional to the 
collision frequency i/m, more collisions increases the friction (Kulsrud and Pierce 1969); 
and b) in a very weakly ionized plasma the collisions are numerous and the whole fluid 
moves together. In this case, the stronger the coUisions the better locking of the gas- 
plasma fluid, and the damping of the wave (which is now proportional to l/t'in) vanishes. 
The Alfven velocity in such a mixture includes the total fluid density m^nj + m„n„. The 
dispersion equation of the Alfven wave is, according to De Pontieu & Haerendel (1998): 



Taking these statements as facts, here we suggest the following item to be added in 
order to complete the physical picture, thus c) the perturbed velocity of the gas-plasma 
mixture may be drastically reduced in a weakly ionized plasma (like the photosphere), 
and, consequently, the wave energy flux becomes very small. 

To describe what happens in reality, we here assume the same magnitude of the mag- 
netic field perturbation as in the ideal case discussed above, i.e., taking it as 1%. Due to 
the perturbed magnetic field, there appears the electric field as described above, and the 

— * 

consequent ion motion in the same direction as the perturbed vector Bi. We denote this 

— * — * 

perturbed E x S-drift velocity of ions by V^. Note that this velocity is the same for elec- 
trons, and that we are speaking about fluid velocities. The relaxation velocity of neutrals 
and ions can be obtained from the following estimates (Mihc 1970). Assume that in the 
starting moment the unit volume of the neutrals have a different velocity Vn. In view of 
the huge difference in mass, we neglect electrons for simplicity. The collision frequency 
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Figure 2: Relaxation of the ion and neutral velocities (normalized to Vc = {PinVn + 
^niVj^/^Pin + Pni)) due to coUisious, for parameters appropriate for the solar photosphere. 
The initial velocities for neutrals and ions are respectively and 1.3 ■ 10^ m/s. 

(see Table 1) is extraordinarily high, of the order of ~ 10^ Hz. Compare this with the 
theoretical gyro- frequency for ions Qi ~ 10^ Hz. Knowing that the wave frequency must 
be much smaller, the frequency ordering that we have here is: 



As a result, we can take the starting/maximal value of the ion velocity Vi and estimate 
for the value it will take within the collisional time. 

In view of the ordering ([6]), the time dependence of the velocities of the two fluids (ions 
and neutral) in relative motion, after the initial movement of the plasma due to electro- 
magnetic perturbations has taken place (regardless of the origin of these perturbations), 
is determined mainly by the friction, and can be obtained from the following equations: 



(6) 




(7) 



Simple combinations of these equations yield two integrals of motion: 




This further yields 



m 



■ exp[-(z/i„ + Uni)t], 
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Table 2: Parameters of waves (wave-lengths A in km and frequencies in Hz) propagating 
through the chromosphere for two different ahitudes h (in km). 



h = 1065 



A 


UJ 


kCa 




0.1 


311 - 1222i 


44855 


3.9 


1 


327 - A5i 


4485 


0.14 


10 


33.1 - 0.45i 


448.5 


0.014 


100 


3.3 - 0.0045i 


44.85 


0.0014 


500 


0.66 - 0.0002i 


8.97 


0.0003 


h = 1990 


A 


UJ 


kCa 


UJi/uJr 


0.1 


69666 - 732i 


69829 


0.01 


1 


6891 - 722i 


6983 


0.1 


10 


371 - 9A.5i 


698.3 


0.25 


100 


36.4 -0.9z 


69.8 


0.025 


500 


7.3 - OMi 


14 


0.005 



— ■ exp[-{Uin + UnM- (9) 

It is seen that the two velocities relax very quickly towards the first term on the right- 
hand sides in these expressions. Taking ~ m„ and = 0, we obtain for the relaxed 
(common) velocity for both species 

Vc = Vi = Vi ~ Vi . (10) 

Hence, if ions start to move due to the electromagnetic force caused by the wave, the 
strong friction will result in a common velocity, which i) is achieved very quickly, and 
ii) is much below the starting velocity of the ion fluid. This behavior is presented in 
Fig. 2 for the photospheric parameters at h = 250 km. 

Here, for the same perturbation of the magnetic field (1%), we have Vc = 1.15-10~^Vi = 
0.15 m/s. Compare this to the velocity in the ideal case Vi = 1.3 ■ 10^ m/s. Note also that 
both Vi and Vc are below/much below the sound velocity = 8.9 ■ 10'^ m/s, respectively. 
Hence, neglecting the pressure (compressibility) effects is justified. Because Vc is so small, 
in Fig. 2 we normalize velocities to Vc and give the plot for f„ and l/fj. It is seen that 
the velocities of both neutrals and ions relax towards the same (normalized) value (= 1) 
within a time interval that is many orders of magnitude shorter than the wave oscillation 
period. As a result, using ([5]), we have the flux in the weakly ionized plasma (for = m.p) 
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miUi + m„n. 



miUi 



'n 



) 



3/2 



(11) 



F ^ 10"^ ■ Fi 



id — 



5.3-10"^ J/{ni^s). 



(12) 



It is seen that the actual flux is always small for any realistic amplitude of perturbations. 
For example, even taking exceptionally strong magnetic field perturbations, e.g. Bi = Bq, 
yields F ~ 5 J/(m^s). Consequently, regardless of the physical mechanism for eventual 
excitation of the Alfven waves in the photosphere, the expected amplitude of the perturbed 
velocity is of the order of 0.1 m/s, and the energy flux of the waves is about one million 
times smaller than the one obtained from the ideal models that assume a perfect coupling 
between the plasma and magnetic field (i.e., ignoring the effects of collisions between ions 
and neutrals, and the consequent weak magnetization of plasma species). 

The estimated flux presented above is obtained for rrii = nip. Taking the more realistic 
value irii = 35mp [TS], we obtain only F = 0.02 J/(m^s). Assuming in addition a stronger 
magnetic perturbation of 10%, we obtain F = 2 J/(m^s) and the common velocity ampli- 
tude f c — 9 m/s. The actual flux may have larger values, e.g., due to stronger magnetic 
field perturbations, but the linear wave theory then becomes unapplicable. 

Since the electromagnetic force still acts on the plasma volume in the time interval 
z/j~l, after the initial movement of the plasma has taken place, one could claim that the 
flux may be higher. Yet, in view of the frequency ordering ([S]), which implies a difference 
of many orders of magnitude, the inclusion of this additional electromagnetic effect in Eq. 
([7]) is insignificant. In fact, it is questionable and indeed unlikely that the ions can really 
achieve the assumed starting perturbed velocity Vi in the first place. This is because the 
assumed value for Vi follows from the ideal case discussed above, with time and spatial 
scales determined by uj~'^ and k~^, respectively, resulting in the characteristic velocity Cq, 
while in the collisional case that we have here, these scales are determined by v^^ 
kJ^, where nj^ = Vn/i^in is the ion mean free path. Hence, the characteristic velocity that 
now appears instead of Ca is Vn = {nTi/miY^'^ and it is about 2 orders of magnitude lower 
than Ca, and a realistic flux should be even smaller than the value obtained earlier. One 
could also argue that the case discussed above, \4 = 0, may look the least favorable for 
the propagation of the wave because neutrals are initially usually in the state of motion. 
Clearly this does not change anything, because in this case, due to the strong collisions. 
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the ions will nearly be in the same state of motion (see in the Sect. 4 below), while 
the ion velocity Vi would still describe an access ion momentum obtained due to the 
electromagnetic perturbation, which neutrals initially do not take part in. 

4 Discussion 

Standard estimates of the wave energy flux through the solar photosphere assume a plasma 
velocity in the photosphere of the order of 1 km/s. This implies two effects: that plasma 
particles move with the observed speed of the convective motion, and that the motion 
of plasma species involves the magnetic field perturbations due to frozen-in magnetic 
field effect. The first effect is only partly satisfied. If in the equilibrium neutrals move 
perpendicular to the magnetic field, say in the x-direction, the plasma particles will move 
also due to the friction effect. The induced velocities of ions and electrons can be calculated 
from Eqs. ([2]) and ([3]) reading 

-t ( ^in -. , ^in \ /I 

Vil.0 = Oi I --^e^ X Vn±0 + Q2^ri±0 I , (13) 



and 



J / ^ 672— # — # ^£2—* — * ^ 671 ^ 6 — # 

Ve±0 = tte I Y^e^ X Vn±0 + — 6^ X Vi±o + -^^nLQ 



-^^m , (14) 



where 



_ 1 



The ion drag velocity (in the x-direction) and the drift component (in the y-direction) 
become, respectively, 

_ _ 1 

i + i ij / 



and 



Vi04rift - Viyo - -TTT— — 2~7o2 - - — ViO,drag- liOj 



The corresponding electron components are 

VeO,drag=VexO = C^eVnx07r- + TT H 1 + 



y2 

in , 



(17) 
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and 




) 



-1' 



(18) 



The induced ion and electron velocities are not necessarily equal, implying the presence 
of equilibrium currents. For the same parameters used in Table 1 and taking the neutral 
velocity of 500 m/s, a.t h = 250 km we have the drag and drift velocities for electrons 315 
and 240 m/s, respectively. The ion drag velocity is almost equal to the neutral velocity. 
This is all due to the fact that the plasma particles are un-magnetized, Qi/ui = 3.6 ■ 10^^, 
Vle/ve = 0.76. However, due to the same reason the frozen-in condition is far from reality 
and the ion/electron motion perpendicular to the magnetic lines does not necessarily 
involve the appropriate movement of the magnetic lines. The actual motions develops as 
described in the previous section. 

In view of the item b) discussed in Sect. 3 [8], such an upwards propagating wave is 
very weakly damped in the photosphere (the damping is proportional to l/uin). This holds 
provided that the wavelengths exceed a certain minimal value. However, it will in fact be 
more strongly damped in the upper layers, e.g., in the chromosphere where the amount 
of neutrals decreases but the damping is proportional to Um- For the chromosphere this 
can be directly demonstrated by solving the dispersion equation that follows from ([Tj) 
where the perpendicular currents are calculated from Eqs. ([2])-(j4|). The expressions are 
very lengthy and we shall not give them here. 

As an example, assuming the wave propagating towards the chromosphere, the dis- 
persion equation is solved for several wavelengths A, with all collision frequencies in- 
cluded, at the ahitude h = 1065 km where T = 6040 K, n„o = 1.71 ■ 10^^ /m^, and 
no = 9.35 • 10^^ /m^, and at the altitude h = 1990 km where T = 7160 K, n„o = 10^^ /m^, 
and no = 3.9-10^^ /m^. The results are given in Table 2. It is seen that shorter wavelengths 
are more damped at lower altitudes. In the same time, longer wavelengths (i.e., those 
that are presumably better transmitted by the photosphere) are in fact more damped at 
higher altitudes. This mode behavior is in agreement with the model of Kulsrud & Pierce 
(1969). However, this trend certainly can not continue because neutrals vanish at still 
higher altitudes. 

We stress that the equilibrium parameters change with the altitude and for the large 
wavelengths the model becomes violated. A numerical approach should give more reliable 
results. Such an approach could help explain how and where the Alfven waves, that were 
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recently detected in the chromosphere [?], are generated. 

Our analysis is based on the presence of a temperature minimum in which most of the 
plasma is neutral, which is predicted by hydrostatic models averaged in space and time, 
such as VAL and FAL. If flux tubes for some reason lack this temperature minimum, the 
analysis we present here may not be an accurate description of how Alfven waves are 
generated in the photosphere. 

5 Conclusions 

The physics of a multi-component weakly ionized plasma, like the one in the solar pho- 
tosphere, is highly complex. Various aspects of this complexity have been pointed out in 
Sects. 2-4. For a temperature of about 0.5 eV, typical for the photosphere, there is a 
plethora of processes that take place and that are nontrivial to include in an analytical 
work like the one presented here. Among others, these include the elastic and inelastic 
collisions, the charge exchange being an important sort of the latter, which imply the 
creation and loss of plasma particles. Yet, in spite of that, some conclusions regarding 
the importance of the electromagnetic Alfven-type perturbations in such weakly ionized 
environments can be made with some certainty. The important conclusion is that if we 
assume Alfven-type waves generated around the temperature minimum, in fact their am- 
plitudes are such that the wave energy flux is very small. The main reason for this is ion 
collisions, which are so frequent that ions almost do not feel the effects of the magnetic 
field. As seen from Fig. 1, in such an environment the ion motion is very similar to the 
Brownian motion of atoms and molecules in a gas. The physics presented here should 
be taken into account in the estimates of the role of the Alfven waves generated in the 
solar photosphere in coronal heating scenarios. However, the solar photosphere is only a 
thin plasma layer and the parameters in the solar atmosphere change with the altitude 
and so does the physics of the Alfven waves. Our analysis suggests that if these waves 
are generated below the chromosphere, they cannot probably be generated around the 
temperature minimum, but perhaps would have to come from lower down, i.e., below the 
surface where the plasma is again much more ionized and the ion-neutral collisions are 
not significant. 
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